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ABSTRACT 
The existence of the two classes of thermal neutrons about the 

Bragg cut-off in a crystalline structure has been experimentally 
verified by the use of dysprosium-aluminum wire detector. The space 
dependent neutron flux in a 24x24x48 inch rectangular parallelopiped 
graphite pile, GLC grade HC and density 1.60 Tyla has been investi- 
gated at 295°K and 481°K. At each temperature a distinct "knuckle" or 
change of neutron distribution is readily apparent, although this 
degree of change decreases with increasing temperature because the 
neutron cross section below the Bragg cut-off increases with increasing 


Cemperature. 
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I. INTRODUCTION 


A distinct characteristic of the neutron scattering cross section 
is known to exist in crystalline materials such as beryllium and 
carbon. The cross section above an energy level corresponding to the 
Bragg cut-off is relatively constant. However, below this energy level 
the cross section experiences a sharp drop, dividing the thermal 
neutrons into what may be called the cold neutron and thermal neutron 
groups. Additionally, within the cold neutron group, the cross section 
varies with the temperature of the material for any specific neutron 
energy. The existence of these groupings is graphically illustrated by 
the variation of diffusion length across the thermal neutron energy 
spectrum. The temperature dependence of the cold neutron group has 
been studied analytically and calculations made to determine the 
diffusion lengths in graphite at temperatures of 300°, 200°, 150°, and 
100°K by Ahmed [2]. 

It was felt that experimental verification of this theoretical 
temperature dependence was in order. A thermal pile was constructed 
with provisions for raising the graphite temperature up to 650°K. 
Although Ahmed's calculations were based on below room temperature 
effects, the same dependency exists for the higher temperatures, as 
can be seen in Table 1 which includes values calculated by Trikha, 


Beeal. {21}. 


TABLE 1 


TEMPERATURE (K°) THERMAL GROUP Ly. COLD GROUP Ly. REFERENCE 
100 2252 lel s> on 
150 ANG, 45.00 2 
200 62.80 - 43.10 2 
300 56.75 24.25 21 
300 58.30 25.00 2 
373 Dons 20.38 aa 
a23 Goi? 15-96 21 
609 SS: 14.49 21 


II. THEORETICAL CONSIDERATIONS 


A. CRYSTALLINE EFFECTS ON NEUTRON SCATTERING 

Since the wavelengths of thermal neutrons are of the same order of 
magnitude (about 1.8 A at room temperature) as the interatomic plane 
Spacing for crystalline materials such as graphite, certain optical 
effects are observed. That is, all the interference effects associated 
with light and x-rays can be duplicated by the neutrons. It will be 
shown that there are energy regions where these optical effects may or 
may not be present. When they are present we have coherent scattering; 
that is, the kind of scattering that produces a wave capable of inter- 
fering with the incident neutron wave. In the absence of optical 
effects, the scattered wave does not interfere with the incident wave, 
resulting in "diffuse" scattering, called incoherent eavescinn, 

The coherent scattering may be composed of both elastic and.. 
inelastic scattering. Incoherent scattering may also be either elastic 
or inelastic, broken down cates day as follows: elastic incoherent may 
be due to the presence of different isotopes in the atomic lattice 
(which gives diffuse scattering), or due to the neutron spin relative 
to the nuclear spin; inelastic incoherent which is due to the energy 
exchange between a neutron and the medium's lattice vibrations, with 
the neutron absorbing one or more phonons (the lattice vibration 
quantum) of energy. This phonon energy may be from vibrations in =z y 
of the degrees of freedom. 

The total cross section may then be written for the thermal region 
as: 


0 = o0 +0 
a 


+o + 
t oe ie 


1 dis 


where Oo. is the capture cross section, oF is the coherent elastic 


scattering cross section, o is the incoherent elastic scattering 


dis 
cross section due to isotope presence or neutron/nucleii spin, and 

oy is the inelastic incoherent and coherent scattering cross section 
due to phonon interchanges. In the thermal region, only the one phonon 
effect (neglecting the two or more phonon effects) makes any contribu- 
tion to inelastic scattering [13]. 

In the case of graphite some of these terms may be eliminated. 
Since coherent inelastic scattering is low compared to the coherent 
elastic scattering it is generally neglected at energies below 1.0 ev. 
The capture cross section of graphite is small but proportional to 1/v. 
Since graphite is generally monoisotopic and has zero spin, the 
incoherent elastic scattering term drops out. Therefore there exists 
only 


So = oo to (above Bragg cut-off) 


or 
oc = ogo, to (below Bragg cut-off) 


In coherent scattering the mechanism by which waves capable of 
interfering with the incident neutron wave is governed by Brageg's law, 


written as 
nA = 2d sin 8 


where n is the order of reflection, } is the wavelength corresponding 
to the neutron energy, d is the plane spacing and 6 is the angle 
between the direction of neutron propogation and the crystal plane. 
From this basic equation it is seen that there are many wavelengths 


(or energies) for which the coherent scattering mechanism is possible, 


based on many different plane spacings. There is also a minimum energy 
corresponding to a maximum wavelength, Aax = 2d, ree which 
coherent scattering is not possible for any given plane spacing. For 
the largest plane spacing, the 002 plane in graphite, this minimum 
energy is called the Bragg cut-off, and no further coherent scattering 
can take place at lower energies. For graphite, the average room 


temperature crystalline 002 plane spacing d is 3.3538 A (Angstroms) . 


Using the deBoglie relation 
— i ee 
P mv 


or h 


vo= Fm 


where h is Planck's constant, p is momentum, m is the neutron mass and 
vy is the neutron velocity, the energy of this Bragg cut-off can be 


determined. 


Substitution of the values for neutron mass, Planck's constant and 
plane spacing yields a Bragg cut-off energy of 0.001823 ev. Thus, 
for all energies above this value (a region occupied by neutrons of 
the thermal group) coherent elastic scattering is the dominant factor, 
while at lower energies, for neutrons of the cold group, incoherent 
inelastic scattering alone dominates. 

Although the frequency of the lattice vibrations is largely 


independent of temperature, the amplitude increases with increasing 


temperature [8]. This increased amplitude results in increased lattice 
energy and therefore larger one phonon energies, giving greater 
incoherent inelastic scattering. Additionally it can be shown [10] 
that the incoherent inelastic scattering is proportional to the wave- 
length, or follows the "1/v law", and that it increases rapidly with 
temperature (about linear by observation, as opposed to the r behavior 
predicted by the usual Debye frequency spectrum). Capture cross 
section, however, is temperature independent for all thermal energies. 
Figures 1 and 2 show the asihonant and incoherent scattering cross 
section regions. The combined effects of these graphite scattering 
cross sections, along with its capture cross section is seen in 


Figure 3. 


B. DIFFUSION OF NEUTRONS IN A CRYSTALLINE MODERATOR 


The diffusion length can be defined by 








2 l ; af 
Ht 7% Nee i 2, 
a 
L2 = i ° LE mae 
3 no »Y 
tr a 


where Nake is the transport mean free path, dS is the macroscopic 
Capture cross section, n is the number of atoms per unit volume and 
ey is the transport cross section. In the thermal neutron group 

region, above the Bragg cut-off, the value of Oy will be largely 
that of the elastic coherent scattering cross section, oe It is 


generally assumed to be given, away from the immediate vicinity of the 


Bragg cut-off, by its asymptotic value O. (1 - =) [1}, where o. 
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is the microscopic (coherent elastic) scattering cross section and A 

is the atomic mass number. The temperature dependence of this cross 
section is very small. In the cold neutron group region O., can be 
taken to be equal to oO, >» the incoherent inelastic cross section due 

to the one phonon effect. As pointed out above, and shown in Figure 2, 
this cross section increases with temperature. Therefore, the diffusion 
length should decrease with increase of temperature for energies in the 
cold neutron group, while remaining essentially constant with tempera- 
ture for energies above the Bragg cut-off. 


The time independent, isotropic diffusion equation for thermal flux 


can be written as 


v9G) -+5 6G) = 0 (1) 


where 


and where D and are the average values in the thermal region for 
diffusion coefficient and macroscopic capture cross section, respectively. 
Equation (1) can be written for two regions of energy, the cold and 

thermal groups, with the values of D and Dm being determined over 


the appropriate energy range. 


Eo ~lev 
cs DCE) $(E) dE 
0.0018 


.0018 
Cay ~lev 


vl 


oCE) dE 
0.0018 
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Thus , 


and 


J D(E) $(E) dE 
AL) | See eee 


2 0.0018 
J $(E) dE 
0 
Elev 
» , F) ¢(E) dE 
z cA 0.0018 
al 
Ey lev 
O(E) dE 
0.0018 
0.0018 
J ~ (E) o(E) dE 
= 0 © 
v2 ane 


0.0018 
J ¢(E) dE 
0 


9) -—54,@ = 0 (1-a) 
Tl 

V $2) ~~ 4) = 0 (1-b) 
T2 


where the subscripts 1 and 2 refer to the thermal neutron group and 


cold neutron group, respectively. With the source placed at the center 


of the base of a rectangular parallelopiped, $¢(x,y,0) will be 


Symmetric about the origin and the flux will be an even function of 
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both x and y throughout the pile. The solution to equations (l-a) and 
(1-b) has the same form. It can be shown that (dropping subscripts for 


simplicity) the solution is [15] 





(x,y,z) = 2 An sinh Yoon 672) : oe) So * cos ara (2) 
mn : 
odd 
where 
a/2 f@a/2 
es 4 mmx nt 
An = $,.(x,y,0)cos—= cos— > dx dy 


ae sinh(y no 
is ~a/2e-a/2 


(a is the relaxation length and is expressed as 


ae = 7 +h? + ot 
T 
c is the vertical extrapolated height of the pile. 
a is the extrapolated horizontal dimension of the pile. 
X,y,Z are the co-ordinates from'an origin at the bottom center of 
the pile. 


Since measurements are made along the centerline of the pile 


where x = y = 0 , equation (2) simplifies further to 


o(0,0,z) = ¢f fies sinh(c-z) 
m,n 
odd 


or in exponential terms 


yy. (e-z) —y__(c~-z) 
$0.2) = > 5 A Te™ ce) (4) 
m,n 
odd 
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In a pile measurement the height c is much larger than the diffusion 
length so that in areas not too near the top, the term Y an 6072) 


becomes very large and the second term of equation (4) can be neglected. 


Then 
Peal Yon® . Vann” 
$(0,0,z) = 5 se Aun e e 
odd 
and by defining 
C = a A 9 mn 
mn 2 m,n m 
odd 
ae (5) 
¢(0,0,z) = = Cae -_— 
odd 
=a 2 Nae & —Y_ 2 
ss mn mn mn 
¢(0,0,z) Cia e + C13 e + C25 e +... 


Hughes [10] points out that as the distance from the source increases, 
the effect of the higher harmonics (Y¥43> Ya. et cetera) becomes 
negligible and that the neutron distribution is practically exponential 
beyond that distance. Also Kaiser and Kimel [12] have computed the 
relative contributions of the ¥43 harmonic for a graphite pile as a 
function of distance from the source. They show 43 to be approxi- 
mately one-tenth the contribution of Yq for a constant thermal 
source at a distance of 20 inches from the source. 

Therefore by concentrating on the mid-axis range — from 35 cm to 
80 cm — it is possible to neglect the higher harmonics and end effects 


so that 
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ooo. 2) = C4 © si 


or 
Ind = In C14 Taine (6) 


which is a simple, one degree polynomial equation of intercept ln Cia 


and slope ona Ing versus z plot. The form of this 


ial 
equation is valid for both thermal and cold neutron groups. 

Least squares polynomial approximations over the region of validity 
for each group results in the determination of their respective Yau1° 
Since the relaxation length is a function of Ly and it has been shown 
that diffusion length for an infinite medium changes across the Bragg 
cut-off, it follcws that the values of Yq, OF the slopes of the plot, 
should change also. This will occur at some point in the pile at which 
the cold neutrons have been filtered through in appreciable number. 

Similarly the degree of slope in each instance will change with a 
temperature change. The greater degree of change will be present in 
the cold group region neutrons, that is, the region where the number 
of cold neutrons is appreciable. Since, as was pointed out, the inco- 
herent inelastic scattering in this region is highly temperature 
dependent, an increase in temperature will result in increased 
scattering cross section and reduced diffusion length. A lower 


diffusion length is associated with an increased slope (refer to 


‘an 
equation (3)). In the thermal neutron group, the transport mean free 
path has been found by Lloyd, et al. [16] to decrease slightly with 


temperature. Therefore with higher temperatures the diffusion length 


would increase and the slope Yu would decrease. The overall effect 
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of high temperature will then be to increase the slope of the first 
part of the neutron distribution curve (before the "knuckle") and 
decrease the slope of the second (or cold neutron) part of the curve 
(beyond the "knuckle"). At a sufficiently high temperature the 
"knuckle" will vanish and a single slope will exist for the entire 
thermal range. According to cross section data from Ref. (Spee chars 


temperature should be higher than 1020°K. 
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III. PROJECT DESCRIPTION 


It was proposed to verify experimentally: first, that the thermal 
neutron diffusion lengths above and below the Bragg cut-off energy did, 
in fact, differ; and second, that this difference decreased as tempera- 
ture increased. The neutron flux plot should change from a distinct 
"knuckle" at the location where the number os cold neutrons becomes 


appreciable to a single, continuous line at very high temperatures. 
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IV. TEST APPARATUS 


A. DESCRIPTION 
1. General 

It was originally planned to utilize the AGN 201 reactor at the 
Naval Postgraduate School reactor facility as the neutron source for 
this experiment. However, the physical Bone reuctdon of the reactor did 
not provide for removal of any segments of its shielding. It would 
have been necessary to construct the graphite pile outside the reactor — 
a distance of some 7.5 feet from the reactor core center — over one of 
the access holes or the "glory hole’. It would be necessary for the 
reactor to be operated at its maximum power of 1000 watts (which could 
only be sustained for about ten minutes) in order for sufficient neutron 
flux to be obtained in this manner. This short time period was quite 
insufficient for any flux measurement in the proposed sigma pile. 

Therefore it was decided to utilize three Plutonium-Beryllium 
neutron sources: two sources of 1.8 x 10° n/sec each, and one source 
of 5.4 x 10° n/sec. These were the only other continuous neutron 
sources available. The most probable energy of these neutron sources 
is about 4.0 Mev [3], so it would be necessary to provide for additional 
slowing down time. 

2. Graphite Pile 

The graphite pile consisted of a rectangular parallelopiped 
constructed of one hundred twenty eight (128) graphite bricks, each one 
3 x 6x 12 inches. The base was two feet square and the height was 
four feet. Figure 4 gives a cross sectional view of the pile's general 
arrangement. Great Lakes Carbon Corporation, Grade HC, graphite of 


density 1.60 alias was utilized. The use of this relatively impure 
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graphite as opposed to AGOT nuclear graphite was considered acceptable 
in view of the experiment's primary aim in investigating temperature 
dependence rather than absolute values. A three sixteenth inch vertical 
hole was drilled down the entire four foot centerline of the pile to 
provide for the detecting wire used. 

3. Heating and Insulation 

Heating of this pile was achieved by twenty calrod strip heaters 
of 21.5 inches heating length, 1.5 inches width and 3/8 inch thickness 
arranged about all four sides, the top and bottom. Close fit grooves, 
1.5 by 3/8 inch, were cut in the graphite bricks, three inches in from 
the side faces and six inches in from the top and bottom, prior to 
assembly to accomodate these heaters. The grooves were a full 24 
inches long, the same size as the pile, in order to provide for thermal 
expansion of the heater elements. Figure 5 shows the external arrange- 
ment of graphite bricks and groove openings. Each heater was rated for 
750 watts at 220 volts, providing for 15 kVA of heating capacity. 
Heaters were regulated in pairs by ten variacs which provided excellent 
control of heat throughout the pile. Figure 6 shows control table 
arrangement with the variacs and thermocouple readout instruments. 

The design of insulation was entirely dictated by the availa- 
bility of insulating materials within the campus. Fire brick was used 
as the base insulation since it was strong enough in compression to 
support the 1590 pound weight of the graphite and also had excellent 
insulating properties. A sheet of three sixteenth inch asbestos cement 
board separated the fire brick from the half inch steel plate which 
acted as a table to support the pile over the neutron source. Both the 


steel plate and fire brick had a three inch hole at the pile centerline 
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Control Table Arrangement 


to reduce neutron scattering in that immediate vicinity. Two types of 
insulating block, using readily available material, were installed 
about the four sides at a distance of about three inches from the pile. 
A two inch layer of glass wool insulation attached to the underside of 
a plywood board encased the top of the insulating walls, also with an 
air gap, thus completely insulating the graphite pile. The air gaps 
were provided in an attempt to reduce to a minimum the backscatter that 
may have occurred from the insulating material, thus making possible 
the use of the simplified diffusion equation given as equation (1). 
Figure 7 shows the lower corner of the graphite pile to illustrate the 


base and side insulation, as well as the heating element connections. 


4. Temperature Determination 


Ten copper-constantan miniature thermocouples were arranged in 
groups of five each in planes that were twelve and thirty six inches 
respectively from the base. In each plane one thermocouple was at the 


centerline of the pile and the other four were three inches from the 
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centerline. The closest any thermocouple was to any heating element 
was 6.72 inches and the most distant was 7.51 inches. Due to symmetry 
of all heating element locations, the thermocouples were equally 
affected by the aggregate heaters. This arrangement of thermocouples 
allowed determination of temperature uniformity throughout the volume 
of graphite for which the diffusion length was measured. A multiple 
switch was utilized to give a common reference junction of 0°C to all 
thermocouple readings. A plain water ice bath in a thermos jug served 
as the reference junction. A digital potentiometer provided the voltage 
readout continuously to within a tenth of a millivolt. A galvanometer- 
type millivolt potentiometer with accuracy to *1.0 microvolt was 
utilized for final voltage determination after the system had settled 


out and for accurate temperature determination during irradiation runs. 
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It was recognized that copper was probably the worst possible 
constituent to use in a nuclear radiation experiment for temperature 
determination. As pointed out by Ross [18], the percentage change of 


resistivity with temperature in a neutron flux can be expressed by: 


Ze” =O OE ot 
ele 


compared to that for constantan: 


hap se Moe ae 


oe 
ee 
II 


where Ap is the change of resistivity, Py is the original resistiv- 
ity, ¢ is the neutron flux and t is the time in the neutron flux in 
seconds. However, the low neutron flux available for this experiment 
and the relatively short exposure times permitted the neglect of this 
shortcoming. As with the decision to use plain water versus pure water 
in the ice bath, the accuracy of termperature determination was deemed 
sufficient to within five degrees centigrade. 
5. Neutron Source 

The three Plutonium-Beryllium (Pu-Be) sources were grouped 
together in a horizontal plane beneath the base of the pile directly 
under the centerline axis. A single one and three quarter inch layer 
of paraffin was placed above the sources to aid in slowing down the 
high energy neutrons before they entered the graphite. Numerous 
arrangements of source height position and thickness of paraffin were 
tried before the above set up was decided upon. It was found that the 
paraffin not only slowed down but absorbed the neutrons so much that 


the flux was reduced to an unacceptable level. Therefore only one 
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layer was utilized. This layer and the source under it were placed as 
close to the base as possible since it was also postulated that the air 
gaps present, and the different medium through which the neutrons 
travelled before reaching the graphite all contributed to the low flux 
within the graphite. Since these neutrons were of a higher energy than 
normally encountered in experiments of this nature, it was thought that 
they would Sg Ae a much larger amount of moderating material to reach 
thermal energies. However a quick calculation of the number of collisions 
required to achieve this energy in hydrogenous moderators (assuming the 


same logarithmic decrement as water) showed as follows: 


in =—2 = 1, 4x 10° 
E 0.92 1 


ea 


4.0 Mev neutrons: 


= 
i 
mt | 


n = 16.5 collisions to reach 1.0 ev 


E 6 
Or ee ae ee * 
2.0 Mev neutrons: n = : ln 5 0.92 ln 1 
n = 15.78 collisions to reach 1.0 ev 


With a transport mean free path of 0.395 cm for paraffin [20] (the 

average distance between collisions), this indicated th=t only an 

additional 0.284 cm would be required to achieve thermal energies. 
Blocks of paraffin surrounded the sides of the source to a 

thickness of two and one half feet and were arranged outside the 

graphite insulation to protect personnel in the area. The neutron 

leakage through this arrangement was surveyed and found negligible. 

6. Neutron Detection 


It was originally planned that a BF, neutron detector, in 


3 


conjunction with a multi-channel analyzer, would be used to take an 


26 


integrated count at each axial position. Investigation showed that this 
detector would withstand temperatures only as high as 150°C. Similar 
restrictions were found to exist for solid state detectors as well as 
indium and cadmium, which melt at 156°C and 320.9°C, respectively. A 
few solid state detectors were available commercially but to achieve the 
sensitivity required for the low flux measurements in this experiment, 
their size was much too large. Gold foils, though commonly used for 
high temperature measurements offered such a long half-life as to be not 
feasible from irradiating time considerations. A detecting wire of 11.4 
per cent dysprosiun-aluminum alloy was finally chosen, based on its high 
melting temperature and convenient (2.35 hour) half-life. A purity of 
99.94 per cent was certified by the manufacturer and impurities listed 
by name indicated no need for their consideration or possible effect on 
the dysprosium capture cross section; although this factor was taken 
into account in computing the total number of dysprosium atoms, described 
in Appendix A. An average atomic weight of 162.5 was determined based 
on the percentage of the various isotopes to be found in natural 
dysprosium. A value of 930.0 barns was used as the capture cross 
section, as recommended [7]. This coincides with the value obtained by 
weighting the cross sections of each isotope in accordance with their 
abundance in natural Dy, Dysprosium has a small response to epithermal 
neutrons [14,19]. The maximum epithermal contribution to the total 
dysprosium activity has been found to be less than two per cent [4]. 
Therefore the activity determined with this dysprosium wire was consid- 
ered to be completely due to thermal flux. The detecting wire was 
positioned in the three sixteenth inch axial hole along the entire pile 
height during irradiation. A thin glass tube held the wire segments in 


position. 
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On initial irradiation tests the reactor facility counting 
equipment was utilized. This consisted of a thallium-activated sodium 
iodide crystal and photo-multiplier, a pulse height analyzer for total 
decay count determination, and a multi-channel analyzer for energy 
spectrum readout and analysis (as a side benefit and not directly con- 
nected with the object of the experiment). It was later found to be 
More practical to utilize the counting equipment of the laboratory adja- 
cent to the experimental room. This involved a Geiger-Mueller tube 
operating with 1400 volt D.C. bias, and pulse height analyzer. The 
efficiency of the G.M. counting system for beta counting was determined 
to be approximately twenty-seven per cent, while oy: of the photo- 
crystal for gamma counting in the reactor facility was about twenty-nine 
per cent, so that no appreciable advantage could be gained by the 


latter's utilization other than familiarization with another piece of 


equipment. 


B. OPERATION 

Test runs were made at room temperature and at an elevated tempera- 
ture of 481°K. The following operational procedures were followed in 
both cases except for temperature control. 

Two main fifty ampere circuit breakers were turned on to provide a 
maximum of 6.8 amperes to each of ten variacs, which in turn provided 
3.4 amperes to each of two heating elements. The variacs were connected 
such that five of them controlled the bottom ten heaters and five of 
them controlled the top ten heaters. Similarly each variac controlled 
the sane Te elements at a particular height so that even ccntrol of 
temperature could be achieved. All variacs were gradually increased to 


Slowly raise the pile temperature to the desired level with the minimum 
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amount of "hunting". It was found that the temperatures at each thermo- 
couple grouping were consistently even, within 2°C of each other. It 
was thus only necessary to adjust the top with the bottom groupings. 
This was a fairly easy task and temperature variations between any of 
the ten thermocouples amounted to a maximum of plus or minus 4.56°C. 

The final test temperature value is the mean of the ten readings, which 
are about constant, over the entire irradiation time span. 

Once the specified temperature equilibrium was achieved, the three 
sources were placed in position under the graphite pile. In order to 
utilize the counting equipment and obtain the activity at a particular 
location it was necessary to cut the detecting vicerinte one inch 
segments. Due to the cost of the wire, segments had to be used over 
again, necessitating a requirement for being able to insert and remove 
segments to a known axial position. This was achieved through the use 
of a thin glass tubing the length of the graphite pile, into which the 
wire segments could be inserted. The bottom of the tube was sealed with 
a glue to prevent segments from dropping out. A maximum of eight seg- 
ments could safely be inserted in the glass tubing without the glass 
fracturing. Since 48 data points were being taken this meant at least 
six irradiation periods for each test run. An average irradiation 
period of eight hours provided ninety per cent of full saturation. 
After irradiation the tube was withdrawn, the bottom two elements 
removed by snapping off the glass tubing at that point, adding a drop 
more glue to the new bottom of the tube, and reinserting the tube, 
maintaining segment axial position, to continue the irradiation proce- 
dure while the two elements removed were being counted for decay 


activity. Due to the low flux, long (thirty minute) counting times per 
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segment were required. Since even during the short time of absence from 
the flux for segment removal a significant decay occurred in the re- 
Maining segments, it was necessary to derive the iterative formula for 
the activity at each subsequent removal time as described in Appendix A. 

Upon removal of the two segments a period of seven minutes was 
allowed to elapse in order to reduce the 75 second half-life of eee 
to a negligible effect. During this time, ten minute background counts 
were obtained, to be averaged together and applied to the segment total 
activity count. During the initial stages of the experiment, using the 
reactor facility equipment, the energy spectrum of the decaying wire 
segment was determined from the multi-channel analyzer. Although the 
decay rates were very small, a typical dysprosium gamma ray spectrum 
was obtained with no indication of any other material being present, 
validating the stated composition of dysprosium and aluminum. ee ache 
out the entire experimental testing, accurate determination and recording 
of time was maintained. 

The individual wire segments were used in the same axial position 


for both runs to insure consistency in obtaining data. 


C. FAILURES AND MALFUNCTIONS 

The scope of the original project definition had to be curtailed 
drastically for a number of reasons, some already stated: 

1. Detector Restrictions 

The high temperatures sought to be achieved drastically 

restricted the choice of detectors available. The state of the art in 
neutron detectors, especially in high temperature ranges, is geared 
toward very high flux ranges as found in reactors. mer recut red the 


use of a dysprosium wire which, while offering negligible self-shadowing 
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and being sensitive mostly to thermal neutrons, also was too small to 
make up the adverse effects of low flux. Since the wire was very 
expensive it was necessary to find some method of re-using the segments 
after they had been previously used. 

A tubing arrangement was decided upon as offering the simplest 
method of positioning and handling the Segments. Glass tubing was 
readily available but, due to its brittleness, there was a high proba- 
bility of tube fracture and subsequent loss inside the pile. This 
method was to be used only as a last resort. Three different metal 
tubes of various sizes were tried. However, in each case the segments 
bound up inside the tubing, necessitating the cutting of the tubes to 
remove the segments. It was finally resolved that the wire was too 
soft to be stacked inside a tube or pushed through it by a rod. The 
slightest buckle produced a jamming against the tube walls, inducing 
further binding. Therefore the glass tubing was resorted to, as 
described in the operating procedures, with utmost caution being taken 
to avoid fracture inside the pile. 

2. Heat Considerations 

The Pu-Be source provides neutrons from the alpha particles 
(given off by plutonium decay) interacting with the beryllium to give 
off neutrons (Be(a,n) reaction). There is the possibility that some 
alpha particles would pick up electrons to become helium gas before an 
interaction with beryllium could occur. Thus it was decided to keep 
the sources at room temperature to prevent a possible rupture of the 
source cladding due to gaseous expansion, despite assurances from the 
manufacturer that no hazard existed. Coincidentally, by maintaining the 


sources at one temperature for all tests a more realistic comparison 
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of data would be possible since the source spectrum could be a function 
of temperature. 

This tremendously complicated the construction of the experi- 
mental set-up. With the existing insulation installed below the base, 
a design was provided to have an air channel over the source and its 
moderating paraffin layer. A small blower motor was to have been 
installed to force cooling air through the channel and out an exhaust 
ducting. However, when it became apparent that the flux level was 
going to be too low, it was necessary to do away with the air channel 
and move the paraffin and source as close to the base as possible. A 
remote reading bourdon temperature gage was installed between paraffin 
layers in order to keep the source temperature under constant surveil- 
lance. 

During heating operations the temperature variation between all 
measuring thermocouples was easily kept within the limits desired. It 
was found however that over a period of time the steel plate began to 
heat up to such a degree as to affect the paraffin. This was traced 
to the three inch square hole in the firebrick below the graphite. 
Although there was also a hole in the steel plate at that location, 
the single sheet of asbestos sheet cement did not offer long term heat 
insulation properties and conducted the heat throughout the steel plate 
underneath. By filling this hole in the firebrick arrangement with 
plastic refractory and a small sheet of insulating block, it was 
possible to eliminate the problem of the hot steel plate. As an 
additional measure, some of the surrounding paraffin in close proximity 
to the steel plate was removed, while leaving the paraffin around the 


source. This provided a cavern effect around the center column of 
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paraffin and source (which remained in place, close to the graphite 
base). The air blower originally planned for was then installed to 
circulate air through this cavern and around the source, rather than 
over it as was originally planned. This proved quite satisfactory and 


no further heat problems were encountered. 


Re. 


V. DATA ANALYSIS AND RESULTS 


A. DATA REDUCTION 

A total of forty eight data points were collected for the room 
temperature test. After the break in the slope was located from this 
run, only the center twenty two positions were investigated for the 
high temperature test. The saturation activity and corresponding 
thermal flux was computed as shown in Appendix A. The uncertainties 
are described in Appendix B. 

The natural log of the flux was then plotted versus the height in 
order to determine where the break in the slope might occur. As noted 
in the theory, only the center portion of the pile was to be investi- 
gated. A general trend or indication of where the fundamental slope 
Yu became predominant was also of interest in choosing points for the 
least squares approximation, however, so all points were initially 
plotted. These are shown in Figure 8. Once the apparent Bragg cut-off 
break point was fixed upon, the points making up the two slopes were 
submitted to separate least squares polynomial approximations to obtain 
a mean value of slope and a standard deviation. By applying Chevenet's 
criterion [9] for rejecting an experimental reading, certain points 
were eliminated for each slope. The remaining data points, those 
within the criterion, were submitted to a final least squares fit. The 
slopes obtained in this manner were entered in equation (3) to determine 


the diffusion lengths. The slopes and diffusion lengths L,, are 


1G T 


shown in Table 2 for both temperature conditions. 
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Fig. 8 
Thermal Flux versus Height for Graphite Temperature of 295°K 
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Fig. 9 
Thermal Flux versus Height for Graphite Temperature of 481°K 
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TABLE 2 


¥4, (relaxation length) 


Temperature Thermal Cold 
a =o al 
Zoo kk 0.077 + 0.013 cm 0.026 + 0.004 cm 
a -1 =~] 
481°K 0.074 + 0.024 cm 0.037 + 0.008 cm 


Ly (diffusion length) 


Temperature Thermal Cold 
295°K 32.72 em Imaginary 
481°K 45.49 cm Imaginary 


B. FINAL RESULTS 

As can be seen from Table 2 the diffusion lengths are much lower 
than those normally found for thermal neutrons in graphite. For 
instance S. K. Trikha, et al. [21] have calculated thermal group 
diffusion lengths of 56.75 cm and 65.17 cm for 300°K and 523°K, 
respectively. Also they have found a cold group diffusion length of 
24.25 cm and 15.96 cm for the same temperatures. In the present case 
the thermal group diffusion lengths are much lower. For the cold 


group, L,, becomes imaginary. That is, the transverse buckling term in 


Ah 
equatican (3) is larger than the relaxation length Yu: DeJuren and 
Swanson [5] also found imaginary diffusion lengths in graphite piles 
with sides less than 28 inches. 


An important, additional consideration is the fact that the Grade 


HC graphite has many more impurities than the AGOT graphite normally 
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used in experiments of this nature. Mantell [17] lists the ash content 
for Grade HC as 0.30 per cent whereas for the AGOT graphite, the ash 
content is 0.07 per cent. Although no exact content breakdown listing 
was found for this particular grade, other manufacturers of graphite 
for similar applications list varying degrees of iron, vanadiun, sulfur, 
calcium, titanium, aluminum and boron as constituents of ash. The 
presence of these elements could have a drastic effect upon reducing 

the diffusion length throughout the pile. 

Appendix B discusses the uncertainty analysis and it is seen that 
very high inaccuracies are present for almost all data points, due to 
the low total count received and the relatively high background count. 
It should also be noted that the steepness of the thermal group slope 


makes it very sensitive to data variations. 
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VI. SUMMARY 


It has been shown by theoretical analysis that the thermal diffusion 
length as illustrated by the relaxation length does, in fact, differ 
above and below the Bragg cut-off. Further, this difference decreases 
as temperature is increased, indicating that a single straight line 
(slope), or diffusion length, is approached throughout the entire thermal 
range at a sufficiently high temperature Ty? where To is greater than 
1020°K. The slopes listed in Table 2 and shown in Figures 8 and 9 
verify this analysis. Although the slopes thus obtained for the neutron 
group above the Bragg cut-off give lower than Peas LTeheraen group 
diffusion lengths, the trend is clearly brought out. 

In the development of equation (2) we have assumed exponential decay 
and that the transverse leakage has no effect on the total flux spectrum. 
However, as shown by Williams [22], for crystalline moderators there is 
a maximum transverse buckling which invalidates this assumption. 
Theoretically, for the case of graphite, values of buckling for sides 
less than 45 inches will lead to a non-exponential decay of the flux, 
and its energy spectrum along the z direction will change. No discrete 
values of the relaxation length would exist in this case. DeJuren and 
Swanson [5] however reported a diffusion length measured in a graphite 
block with 28 inch sides. Only when their sizes were reduced to 20 
inches did their diffusion lengths become imaginary. As DeJuren and 
Swanson showed in their experiments, the transverse leakage of neutrons 
increases along the z direction due to the fact that cold neutrons with 
longer mean free paths escape easier and therefore spectral equilibrium 


is not maintained. 
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This effect is apparent in the cold group region of figures 8 and 9. 
Due to the small size of the present pile, the spatial decay is no 
longer exponential and equation (2) can no longer be valid. 

The heating capacity of this graphite pile has not been approached. 
Further work in this area at even higher temperatures should be done to 
verify this conclusion, which is based on only one high temperature 
test. It would be recommended however that provisions for a stronger 
neutron flux within the graphite be provided, in which case the thermal 
neutron sensitive dysprosium could still be used. The possibility 
should also be investigated of using a water cooled BF, detector in 
conjunction with a multi-channel analyzer. Either case will reduce the 
uncertainties presented here. Another important point is that the 
block size must be increased beyond the border line between the imagi- 
nary and real diffusion length regions. In fact, the size must be large 


enough to insure the existence of a diffusion length for the cold 


neutron group. 
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APPENDIX A FORMULAE FOR DATA REDUCTION 


In determining the neutron flux, use is made of the radioactive 
atoms of dysprosium which are produced by neutron absorption within 
the wire detector. The rate of production R less the rate of decay A, 
called the activity, is equal to the rate of change of these radio- 
active atoms 


dn 


at R- An 


where n is the number of radioactive atoms and A is the decay constant. 
For steady flux cases R is a constant and the solution to the above 
equation can be written 
-At 
A = jn = AG -e ) (7) 
where A.» the saturation activity, is the equivalent of the rate of 


production R for very long irradiation times. After removal from 


irradiation the rate of decay will decrease exponentially as 


where t; is the time after removal and ae is the activity at removal 

time. If counting is employed to determine the number of radioactive 
atoms n existing over a given period of time, those counts C, less any 
background counts B, will be equal to the difference between the rate 


of decay at the beginning ts and end t, of counting time, divided by 


b 
the decay constant and multiplied by the efficiency e of the counting 


equipment. 
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A A 
| poe 
n= C-B = e€[—+e ra - be beg 


NGHeEEeeeeres [ @ cate. arte 5 (8) 


In the case where the detector wire is alternately removed and 
inserted in a neutron flux, a recursive formula was derived for A,. 
Referring to Figure 10, the initial activity at t = 1 will be simply 


Ay SUA, dls BE sal (9) 


Between t = 1 and t= 2 , the value of A, will decrease an amount 


1 
A, Eo? . Thus the value of X will equal 


X= A, pat oe (10) 
A 
ier Ay 1 3 
i y l 
1 | 
X ene eee ee Ne ns i- > i 
t { 
* tt 1 4 
= t i { ( i 
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Fig. 10 


Repeated Exposure and Withdrawal From Flux 
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Considering an exponential rise again at point X and time starting over 


is 


at t = 2, the value of Ay 


> 
T 


K+ (A -X)C1- e "12 ) 
Nese Ga A. C(1- ap i2 ) (11) 


Combining equations (10) and (11) 


“Nt 


= “A (togtt ja) - 
A, = A, e 1 == A. (l-e ) 
and from equation (9) 
A, = A.C 1-e tiny e*to2tti2) 4 4 (1 - eAti2 y 


Which can be written 


A, = A T1- e424 (1 - tar ye o2*ti2) | 


A similar line of reasoning can be applied to A, 


Ca Be tio “)(to3ttz3) 
A. A. [l-e + A, e J 


the recursive formula being 


A, = AL 1- PS pee a Sa ERY 9 2 ee an 


Equation (9) is still valid for j = 1. From equations (8) and (12) 


ACC - B) 
ent pe ise; Ayy 7) tosttis) ek paca) ob ») (13) 
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For the case where the flux is constant the saturated activity can 


be written 


A = N, fio@om aE 


Ss 


where Ny. is the total number of target neutrons and 0) is the 
capture cross section at energy E. Since the dysprosium was within two 
per cent of indicating only thermal energy neutrons, the integration 
can be carried out over the thermal energy range. Considering a 
Maxwellian flux and assuming a 1/v behavior for dysprosium, the 
saturation activity can be written 


TL 


2 _ oy) 1/2 
A, = Npo¢C Gr) 


ks) 


where T is 293.16°K, T is the medium temperature, and oe 920 b. 
is the capture cross section recommended in Ref. 7. It is based on the 
weighted cross sections for Dy. isotopes by their abundance in natural 


dysprosium. Equations (13) and (14) now give, for the thermal neutron 


flux 
A(C = B 
i (15) 
TT 
=) eee eed S88 Gineely 2 

ak Ce e ) oO. ( iG. ) N,, 

Np» the total number of target neutrons, can be found by 
N. = V(¢ ee) ) v 

. A 


where V is the wire volume, WT is the wire weight, A. is Avogadro's 
number, A is the average atomic mass, and v is the fraction by weight 


of dysprosium in the Dy-Al wire. 
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APPENDIX B UNCERTAINTY ANALYSIS 


As stated earlier the low thermal flux resulted in a relatively low 
count rate compared to the background count. An analysis was conducted 
on each of the data points to determine the accuracy of the resulting 
el tax . 

Considering a total of C counts as the mean value, its standard 
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deviation can be written as C » and the total count would be reported 


as Ct ge » The same criterion is applied to the background count B. 
In obtaining the difference, C - B, over the same time interval, the 


value with its standard deviation is 


= 5) + (0, + 05°) ai 
where the standard deviations are 
o. = rm and 
oO, = pi/2 


This value is converted into a percentage error by 


Zz Zo 1/2 
(9, + OD ) : 
(C-B) + —-———~— 100% 


(c - B) 
Since this was the major source of uncertainty in determining the flux, 
the same percentage error applies to the thermal flux obtained for each 
data point. These uncertainties in terms of percentage are tabulated 
in Tables 3 and 4 at the end of the Appendix. 
When the data points were plotted there were some obviously erratic 


points. It seemed appropriate to disregard some of the more erroneous 
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ones. Therefore a least squares approximation, based on equation (6), 
was obtained on the computer. The standard deviation of the total of 
all data points was also obtained from the best fit straight line. 
Using Chevenet's criterion for rejection of outlying data, a ratio 
d 
(*) of data deviation from its mean value to the standard deviation 
was obtained for each position. Those omnes giving ratios outside the 
criterion were rejected and those less than the criterion were resubmit- 
ted for a second least squares approximation. It is this second 
approximation for slope that is used for the final results and discussion. 
The ratios obtained as above are also listed in Tables 3 and 4. 

The high temperature uncertainty was determined as follows, utilizing 
the voltage readings in all calculations. When the final value was 


obtained, the ccnversion to degrees centigrade was applied. At each 


of the levels, upper and lower, 


_— ab 
Vv = ae 2 V- 
N i=1 1 
N 
- 1 —_ 2 1/2 
eae Nee a 
i=l 


where UG, are each of the five thermocouple readings at one recording 
time and N is equal to the number of thermocouples. Then the mean and 


deviation for the pile at each recording instance is found 


ao 1 a = 

= — (V V 
Vy 5 ( cant ? 

Aa ea 
Oo, = 5 (o, + ch ) 


where the subscripts u and & refer to upper and lower levels, 
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respectively. The overall temperature value and its standard deviation 


is then obtained from the mean values of each recorded run: 


N 
Vn 7 7 “si 
i=l 
N 
o = [ o 2) 12 
m : Ve 
=] i 


Using the above procedure a high temperature of 207.88 + 4.56°C was 


obtained, 
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TABLE 3 
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Temperature: 295°K 
oF 
Height above Flux Error — Chevenet's 
base (cm) (n/cm2sec) (per cent) : Criterion 
35.21 2354.8 8.60 0.05 L773 
37.65 1926.5 11.61 0.5/ NS) 
40.11 1686.2 le? 5 0.14 1e%3 
42.57 311.7 17.30 azo La73 
45.02 152343 12.62 Grol 17s 
47.53 849.6 Za) 2o9 9 73 
49.99 896.9 yA OW as) 7a S18: 2213 
52.43 976.9 21.30 0.14 2.13 
54.85 761.4 22 4.12 Ze3 
D7 .29 798.4 24.95 ic 2.13 
bo. 15 692.7 2ieo> 1.64 2.13 
62.21 J92.9 23.50 0.69 Zaks 
64.63 38 ee, 2020 3.19 Zel3 
67.03 523.9 32,00 Sez Zug 
69.47 492.4 30.90 3.69 2.455 
ea. OT 748.2 24.30 3.82 2an3 
74.35 436.3 34.70 9.26 2503 | 
76.81 563.1 315,50 4.01 Zo 
72.25 439.8 39620 4.59 2.13 
81.67 389.9 45.90 0.64 23 


TABLE 4 


Temperature: 481°K 
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Height above Flux Error “i Chevenet's 
base (cm) (n/cm2sec) (per cent) = Criterion 
32.75 DI1LO.9 10.35 Or52 1.90 
Bo. 21 8612.5 esi 5 s)3, 4.93 1.90 
37.65 2204.9 20.23 3.48 1.90 
50.11 1713.0 2047 3.87 1.90 
42.57 2169.6 W224/ 0.78 Ish, 
45.02 2502.8 15.74 1.66 1.90 
47.53 1671.0 19.00 Og Lon elt) 
49.99 1547.5 24.40 1.04 1.90 
52.43 1204.2 30.50 1.76 2.10 
54.85 5291.3 31.55 3.01 2.10 
By. 29 890.6 30.90 7.79 74148, 
59.75 1055.6 26.50 3.59 2.10 
B22) 1050.0 39.7) O79 2.40 
64.63 802.9 36.30 0.32 2.10 
67.03 740.4 30.45 ee 2.10 
69.47 1022.5 32.10 0.74 74 MKG, 
peed 828.8 Ways 518, 3.52 2.10 
74.35 646.8 115.00 Ja S32, 2.10 
76.81 27 Dail 149.50 1.43 Zee 
Fg eae 466.9 54.90 6.50 2.10 
81.67 (neg) 74.170 Leg dS) 2.10 
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